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Abstract 

X 

The 6 — > s £ + £ process is studied in the minimal supergravity model in detail. Tak- 
ing account of the long distance contributions from the cc resonances, we calculate 
the branching ratio and the lepton forward-backward asymmetry in this model. We 
find that there is a strong correlation between the branching ratios of b — > sj and 
b — > s£ + £~ processes and that the interference effect can change the b — > s£ + £~ 
branching ratio in the off-resonance regions by up to ±15% depending on the rel- 
ative phase between the long and short distance contributions. Using various phe- 
nomenological constraints including the branching ratio of b — > sj, we show that 
there are regions in the parameter space where the branching ratio of b — > s£ + £~ 
is enhanced by about 50% compared to the SM. We also show that the branching 
ratio of b — > s vv is reduced at most by 10% from the SM prediction. 



I Introduction 



Although the Standard Model (SM) of the elementary particle physics is successful in 
explaining almost all experimental results, it is possible that physics beyond the SM 
exists just above the presently available energy scale. Since new physics may affect 
various processes at low energy such as the flavor changing neutral current (FCNC) 
processes of K mesons and B mesons, new physics searches in these processes are 
as important as direct particle searches at collider experiments. A prime example is 
the b — ► s 7 process. Experimentally the inclusive branching ratio is determined as 
B(6 -> 57) = (2.32 ± 0.57 ± 0.35) x 10~ 4 at the CLEO experiment 0. It is known 
that this process puts very strong constraints on various new physics beyond the 
SM, for example two Higgs doublet model and supersymmetric (SUSY) extension of 
the SM. Along with the b — > s 7 process, another important rare b decay process is 
the b — > sll decay. Although only upper bounds on branching ratios are given by 
experiments for various exclusive modes , this process is expected to be observed 
in the near future at B factories as well as at hadron machines. 

In this paper we investigate the b — > sii decay in the minimal supersymmetric 
standard model (MSSM), especially in the minimal supergravity (SUGRA) model. 
The MSSM is now considered to be the most promising candidate beyond the SM. 
In the MSSM, SUSY partners such as squarks, sleptons, higgsinos and gauginos 
can contribute to FCNC processes through loop diagrams. In order to evaluate 
their contributions quantitatively it is necessary to specify how soft SUSY breaking 
terms are generated. In particular, the soft SUSY breaking terms in the squark 
sector become new sources of flavor mixing, and the K°-K° mixing becomes too 
large if the squark mixing is 0(1) and masses of SUSY partners are in below- TeV 
region ||. In the minimal SUGRA model it is assumed that the soft SUSY breaking 
terms are universal at the Planck or GUT scale. Flavor mixing at the electroweak 
scale can be determined by solving the relevant renormalization group equations 
(RGEs) from the Planck to the low energy scale. It is shown that the constraint 

n — 

from the K -K mixing is easily satisfied in this framework since masses of the first- 
two-generation squarks with the same quantum numbers remain highly degenerate 
at the low energy ||. The FCNC processes involving the third generation quarks 
and squarks can receive sizable SUSY contributions due to the large top Yukawa 
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coupling constant. In particular the 6-^37 process has been intensively studied 
both in low energy SUSY Standard Models and in the minimal SUGRA model H-fl . 
It was observed that the SUSY loop effects can interfere with the SM amplitude 
constructively or destructively depending on the parameters on the model whereas 
the charged Higgs contribution is always constructive. For the minimal SUGRA 
model the B°-B mixing and CP violating parameter in the K°-K mixing, ex, has 
been also investigated in Ref. ]/]]. In this paper we consider the b — > s££ process 
in addition to the 6^57 process to see possible implication on the model by 
future experiments. We observe that the predicted branching ratio of the b —>■ 57 
process and that of the b — > s £ + £~ process are strongly correlated and thus their 
measurements are useful to distinguish the SUGRA model from other extensions of 
the SM. 

The b — > s £ + £~ process in the SUGRA model was analyzed by Bertolini et al. 
in Ref. M. Recently this process was reconsidered taking account of the measured 
branching ratio of b —>■ s 7 and the top quark mass in the context of the low-energy 
SUSY models as well as the minimal SUGRA model H,||. In Ref. H it was noted 
that the b —>■ s £ + £~ process is capable to resolve two-fold ambiguity which cannot 
be distinguished from the branching ratio of b —>■ sj. In Ref. |§ a more detailed 
analysis has been done in the minimal SUGRA model. Compared to them, our 
calculation is improved in several points such as (1) numerically solving RGEs with 
whole Yukawa matrices and soft SUSY breaking parameters with the flavor mixing, 
(2) taking account of one-loop corrections in the Higgs effective potential [|l(J in order 
to determine the proper vacuum expectation value through the radiative electroweak 
symmetry breaking scenario [fElH , and (3) including the interference effect with the 
long distance contribution in calculating the lepton invariant mass spectrum of the 
b — > s£ + £~ branching ratio. It turns out that the third effect can change the 
branching ratio in the off-resonance region by ~ ±15% depending on the relative 
phase between the long and short distance contributions. Taking account of various 
phenomenological constraints including the branching ratio of b — > s 7, we show that 
there are regions in the parameter space where the branching ratio of b — > s £ + £~ is 
enhanced by about 50% compared to the SM. We also calculate the branching ratio 
of the b — > s v V process in the minimal SUGRA model and show that the branching 
ratio is reduced at most by 10% from the SM value. 
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This paper is organized as follows. In section [TI|, we introduce the minimal 
SUGRA model and explain new sources of flavor changing in this model. In section 
PTT] , formulas for b — > s£ + £~ decay are given including SUSY contributions. In 
section |IV|, we present numerical results of our analysis. In section [V|, b — > s v V 
decay is discussed. Section [\T| gives conclusions and discussions. Various formulas 
are summarized in Appendices. 



II Minimal SUGRA model 

The MSSM Lagrangian is specified by the superpotential and the soft SUSY breaking 
terms. The superpotential of the MSSM is given by 

W F = e a(3 ' !,■,,()•; U-lr, + fjxjHfCgDj + .//,; ; //i7,7; / - nH?H$\ , (2.1) 

where the chiral superfields Q,U, D, L, E, Hi and H 2 transforms under SU(3)c x 
SU(2)l x U(1)y as following representations: 

Q? = (3,2,i), Z7 i = A = (3,1 A (2.2) 

o 3 3 

L? = (l,2,~), ^ = (1,1,1), 

fl? = (l,2,-i), fl? = (l,2,±). 

The suffices a,/3 = 1,2 are SU(2) indices and i,j = 1,2,3 are generation indices. 
e a f3 is the anti-symmetric tensor with e±2 = 1. A general form of the soft SUSY 
breaking terms is given by 

- £ soft = (m 2 Q )ij q[i q Lj + (mgr)y fijy + (m^)y d* Ri d Rj (2.3) 

4 hj + (m 2 E )ij e* Ri e Rj 
+A 2 X h\h x + A 2 2 h\h 2 + e a/3 (B /jh^h^ + H.c.) 
+e a/3 «Jy + A Dij h a x tf u d* Rj + A Eij h* F Ll e% + H.c.) 

+{-m § B B + \m w W W + \m 6 G G + H.c), 
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where qu, u* Ri , d* Ri , £u, S Ri , hi and h 2 are scalar components of chiral superfields 
Qi, Ui, Di, Li, Ei, Hi and H 2 respectively, and B, W and G are U(1)y, SU(2)l and 
SU(3)c gauge fermions. 

In the minimal SUGRA model the soft SUSY breaking terms are assumed to 
take the following universal structures at the GUT scale. 

( m Q)ij = ( m u)ij = ( m D)ij = m l S ij, ( 2 - 4 ) 

( m l)ij = ( m E)ij = m l s ij, 
A 2 = A 2 , = ml, 

Auij = fuij A x , A Di j = f Di j A x , A Ei j = f Ei j A x , 

m B = m w = m G = M gx- 

With this initial condition, soft SUSY breaking parameters at the electroweak scale 
are calculated by solving the RGEs of the MSSM flT2 |. We first solve the one-loop 



RGEs for the gauge coupling constants taking a>i(mz) as the input and determine 
the GUT scale, M GU t, where the gauge couplings meet. The Yukawa coupling 
constants at M GUT are also calculated by solving the RGEs from m z to M G ut- The 
values of the Yukawa coupling constants at the electroweak scale are obtained by 
taking the quark/lepton masses, the Cabbibo-Kobayashi-Maskawa (CKM) matrix 
elements and the tan j3 = (h®) / (h®) as input parameters. Then we solve the RGEs 
for all MSSM parameters downward with the GUT scale boundary conditions |(2.4 



for each set of the universal soft SUSY breaking parameters (m , A x , M gX ). We 
include all generation mixings in the RGEs for both Yukawa coupling constants 
and the soft SUSY breaking parameters. Next, we evaluate the Higgs potential 
at the electroweak scale and require that the minimum of the potential gives a 
correct vacuum expectation values of the neutral Higgs fields as (h®) = v cos (3 and 
(h®) = wsin/3 where v = 174 GeV. This is known as the radiative electroweak 
symmetry breaking scenario fllTf . The effective potential of neutral Higgs fields at 
the electroweak scale is given by 



V{h\, hi) = (/i 2 + A() \til\ 2 + (/i 2 + A 2 ) \h° 2 \ 2 + {B n h\ h\ + H.c.) (2.5) 



4 



,2 i „/2 



+^^(i^r-i^i 2 ) 2 +Moo P , 

where Vi oop is the one-loop correction induced by the third generation fermions and 
sfermions [|1IJ. The requirement of the radiative electroweak symmetry breaking 



determines the magnitude of the SUSY Higgs mass parameter \i and the soft SUSY 
breaking parameter B. The explicit forms of Vi oop and the condition of the radiative 



electroweak symmetry breaking used in the present analysis are given in [Appendix A . 
At this stage, all MSSM parameters at the electroweak scale are determined as 
functions of the input parameters (tan/3, m , Ax, M g x, sign(/z)). 

With use of the low energy SUSY parameters determined by the procedure 
described above, we can calculate all the SUSY particle masses and the mixing 
parameters. The 6x6 mass matrix of up-type squark is written by 

-C = (ul,u* R )M? (* L ), (2.6) 




1 2 



[m 



LLjij - (Ml 7 M u ) ij + (mj)ij + m 2 z cos(2/3)(- - -sin 2 6 w )5 ij} (2.7) 



2 3 



{m RR h = (M u Ml) i3 + (ml) ij + m 2 z cos(2p)(^sm 2 9 w )5 l3 , (2.8) 

( m LR)ij = ( m RL)tj = - ft cot /3(Ml)ij + (A^ijvsinp, (2.9) 

where M v is the up- type quark mass matrix, i.e. M Vi3 = fujiV sin /3, and u L is the 
up-type component of SU(2) doublet q. In this weak eigenstate, the mass matrix 
is not diagonal at the electroweak scale. The physical mass eigenstate is given by 
diagonalizing the mass matrix, 



uj = (UuYj " l , (2.10) 




U v (Ml) Ul = diagonal, (2.11) 

where ui is the mass eigenstate. The unitary matrix Ujj induces new flavor mixing 
in the up-type squark sector. In a similar manner, we define the mixing matrices 
Ud, Ui for the down-type squark and the slepton sectors. 
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Ill b^s£ + £~ decay 

In this section, we describe the calculation of the branching ratio and the lepton 
forward-backward asymmetry for the b —>■ s£ + £~ (£ = e,fi) process in the minimal 
SUGRA model. We first introduce the effective Hamiltonian which is relevant for 
the b — > s£ + £~~ process O], 



AC 10 

ffc ff ~E Q(Q)Om, (3-1) 

i= i 

where Q is the renormalization point. For the calculation of the branching ratio, 
the following three operators are important, 

°7 = ^— 2 m h {s La a^b Ra )F^\ (3.2) 
9 = ^(s La ^b La )(£Y£), (3.3) 
do = ^(s L a^b LQ )(£n 5 £). (3.4) 



The explicit forms of all the effective operators Oi(Q) are given in [Appendix B 
Throughout this paper we neglect the strange quark mass. 

The coefficients Ci(m w )-Ci (m w ) are determined by matching the full theory 
with the effective theory at the renormalization point Q = mw The coefficients 
Ci(m w )-C 6 (m w ) are given by 

C 2 {m w ) = -A t , Ci(m w ) = 0, (i = 1,3-6), (3.5) 

where At = VtbV t *. Note that there is no SUSY contribution to these values at the 
tree level. The coefficients CV(mw)-Cio(rnw) are generated by one loop diagrams. 
In order to determine these coefficients at the mw scale, we need to calculate photon 
penguin, Z penguin and box diagrams taking account of new contributions in addi- 
tion to the SM diagrams. There are four classes of new contributions in the SUSY 
model; charged Higgs boson(iJ~)-up-type quark loop, chargino (x~)~up-type squark 
loop, gluino (G) -down- type squark loop, and neutralino (x°)-down-type squark 
loop. Cf{mw) is obtained by calculating the photon penguin diagram Fig. |l|(a). 
C$(mw) is also obtained by calculating the gluon penguin diagram Fig. |l](b). There 
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are three types of diagrams which contribute to Cg(mw)', the photon penguin di- 
agram Fig. |l|(c), the Z penguin diagram Fig. 0(c) and the box diagram Fig. 0(d). 
Since we neglect the lepton mass, there is no charged Higgs contribution to the box 
diagram. Cio(mw) is induced by the Z penguin diagram Fig. [l|(c) and the box dia- 
gram Fig. |l](d). The explicit form of each contribution is given in [Appendix C| with 



use of the mixing matrices at each vertex defined in |Appendix D| and Inami-Lim 



functions given in [Appendix E . 

In order to calculate the b — > s£ + £~ decay amplitude, we need the effec- 
tive Hamiltonian at the rrib scale. By solving the RGEs in the leading logarithmic 
approximation (LLA) of QCD, Ci(rrib) can be related with Ci(mw) as given in 
Appendix Fj [|14}j . With this effective Hamiltonian at the rrib scale, we can calculate 
various physical observables. Since the bottom quark is much heavier than the QCD 
energy scale, we can calculate the inclusive decay width as a free bottom quark de- 
cay. This procedure is justified as a leading order approximation of the heavy quark 
expansion ||15||. The b — > s£ + £r branching ratio is given by 



dB(b 



:£ + t 



d s 



c e v) 



a 

4^ 



A* 



V, 



cb 



1 



f ph (m c /m b ) 



wis) 



1 



4mj 



(3.6) 



\C. 



9 + F(s)| 2 ai(s,m s ,m £ ) + |C 10 | 2 a 2 (s, rh s , rh t ) 



4 

+-\C 7 \ 2 a 3 (s,m s ,m e ) + 12a 4 (s,m s ,rfii)Re(Cj(C 9 + Y(s))) 
s 

where s = (p + + p_) 2 /m 2 and p+(p-) is the four momentum of £ + (£~). Here we 
normalize the branching ratio by the semileptonic branching ratio B(6 — > c e u) 
in order to cancel the mjj factor in the differential width. The function f p h{x) = 
1 — 8x 2 +8x 6 — x s — 24x 4 hix is the phase space factor of the semileptonic decay width. 
Kinematical functions oti - and w(s) are given in Appendix Q . As mentioned 
before we neglect the strange quark mass in the numerical analysis. The lepton 
mass is, however, kept since lepton mass corrections are important in the lower end 
of the lepton invariant mass spectrum. By calculating the matrix elements of four 
quark operators 0\-0 & at the one loop level, we obtain Y(s) 



Y(s) 



g{^, s)(3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) 



(3.7) 
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--2(M)(4C3 + 4C4 + 3C 5 



C 6 )--2(0,§)(C 3 + 3C 4 ) + Y; C 



A\ nz 2 + A + ie£ 

1 1 27 1 1 ■ 



9* 



42 2 



(2 + 



4z 2 
s 



In 



4z 2 



1-1/1 



4z 2 



?7T 



for s > 4z 2 



g 111 6 T 27 T gs 



(3. 



-4/4*-l(2+4£) arctan 




for s < Az 2 . 



In addition to these short distance contributions, there are long distance contri- 
butions from the cc resonances, b — > s J/^ — > and 6 — > — ► Although 
we can avoid large contributions from these resonances by cutting the resonance 
regions of the lepton-invariant-mass spectrum, there can be sizable effects from in- 
terference between the short distance contribution and the tail of the resonances. 
The resonance effects in the h —>■ s transition have been investigated in connection to 
the long distance contribution of the b — > s £ + £~ ||T^|-0] as well as the b 



s 7 



191 



process. In Eq. (3.7) following Ref. [16, 17 we have introduced the resonance term 



Y 

± rf . 



3tt 



MiTji -> e + e-)/m 2 b 

^ S - Mf /ml + iMiYi/ml 



(3.9) 



where k parametrizes the b-s-J/ip and b-s-ip' couplings. Its absolute value is de- 
termined from T(b — > J/ipX) and is given by \k\ ~ 1 Hl6, 17|. In general k can have 
a non zero phase. In the following in order to see the effect of the phase, we show 
results with k = ±lQ. In the actual evaluation of the branching ratio the charm 
mass in Eq. |(3.7)| is taken to be the D meson mass ||16|| . The choice of the charm 
mass is not important here since the branching ratio depends on it very weakly. 

Another observable which is expected to be measured with reasonable accuracy 



in future experiments is the forward-backward asymmetry of the lepton fl7 |. In the 



*In principle k can be different for J/tjj and ip 1 . But for simplicity, we have take the same value 
in Eq. |(3.9) . From the experimental data at least we can show that the absolute value is almost 
the same both for J/tp and ip' ■ 
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center of mass frame of the lepton pair, this is defined as 

Uicose) 1( Ws (b - sin-) - j d(cos6) 1( JTs (b -+ sin-) 

Ayb{s) = \ jj (3.10) 

(d(cos9) ^f^(& - sin~) + J rf(cos^) - *£+£-) 

3w(s)^/l -^Cio [s(C 9 + ReY(s)) + 2C 7 ] 
"|C 9 + F(s)| 2 «i + C? a 2 + §C 7 2 a 3 + 12a 4 C 7 (C 9 + Re Y(s))] ' 

where 9 is the angle between the momentum of the bottom quark and that of the 
l + . 

IV Numerical results 

As explained in Sec. ||, the MSSM parameters are determined by solving the RGEs. 
In the minimal SUGRA model, there are five SUSY parameters; the universal scalar 
mass mo, the gaugino mass M g x, the universal A-parameter Ax, the SUSY invariant 
Higgs mass /i, the mixing parameter of Higgs boson B. Using the condition that 
the electroweak symmetry is properly broken to give the correct Z boson mass, the 
theory contains four free parameters, tan/5, m , M g x and Ax as well as the sign of 
/i. We scan the parameters m , M gX and A x in the range of m < 2 TeV, M gX < 2 
TeV and \Ax\ < 5 m for each fixed value of tan/3. We also impose the following 
phenomenological constraints. 

1. b — > sj branching ratio: 

The branching ratio of the b — > s 7 process is given by 

2 

B(b^ceV)\C 7 (Q)\ 2 , (4.1) 

Most important theoretical ambiguity comes from the choice of renormaliza- 
tion scale Q. The branching ratio changes by about ±30% as the scale Q is 
varied in the range of m&/2 < Q < 2rrib- We fix Q = m;, in this analysis 
and discuss the ambiguity associated to the QCD correction later. From the 
measurement by CLEO Q, the inclusive branch ratio is given by, 

1 x 10~ 5 < B(6 -> 57) < 4.2 x 10" 5 . (4.2) 



B(6 



6a 

si) = — 

TT 
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2. From the recent experiment at LEP 1.5 [2IJ, we impose that all the charged 
SUSY particles are heavier than 65 GeV. 



3. All sneutrino masses are larger than 41 GeV . 



4. The gluino mass is larger than 100 GeV. The lower bound of the experimental 
gluino mass is given by Fermilab TEVATRON collider P^| . Since it depends 
on various SUSY parameters, we take 100 GeV as a conservative lower bound Q. 



5. From the neutralino search at LEP f23 |, we impose 

T(Z^xx) <8.4MeV, (4.3) 

B(Z - XX'), B(Z - X'X') < 2 x lO" 5 , (4.4) 
where x is the lightest neutralino and x' denotes other neutralino. 

6. The lightest SUSY particle is neutral. 



7. The condition for not having a charge or color symmetry breaking vacuum [24] . 



Throughout this paper we fix the top quark mass as 175 GeV, the bottom quark 
pole mass as 4.62 GeV and a s {m z ) = 0.116. In Fig. |^, we show C 7 (m b ), C 9 (m b ) and 
Cio(mb), each of which is normalized to its SM value, for tan/3 = 3,30. In these 
figures we do not include the b — > constraint. We can see that in Fig. ^| CV(m&) 
can be quite different from the SM value and even the opposite sign is allowed for 
tan/3 = 30. On the other hand, Cg(rrib) and Cio(rrib) differ from the SM values by 
at most 5% in the whole parameter space for both tan (3=3, 30. In the calculation 
of CV(m&), there is a one-loop diagram with internal stop and chargino which gives 
a large contribution when tan/3 becomes large ||. When chargino has a sizable 
higgsino component, this diagram is proportional to the product of the top and 
bottom Yukawa coupling constants i.e. -^f^a which grows as tan/? when tan/? is 
large. On the other hand, there are no such terms in the calculation of Cg(m&) and 
Cio( m f>)- In fact, the corresponding stop-higgsino diagram in C 9 {m b ) and Ci (m b ) 
is proportional to the square of the top Yukawa coupling constants, namely ™^ -e , 



tWith use of the GUT relation of the gaugino masses Eq. (2.4) and the LEP 1/1.5 constraints 
on charginos and neutralinos (above 2 and 5), a gluino lighter than about 150 GeV is excluded for 
tan (3 > 2. Therefore the precise value of the imposed gluino mass bound is not very important. 
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which does not grow for large tan/3. Indeed C$(mb) and C\a{mb) could be large if 
tan/3 < 1, but within the framework of the minimal SUGRA model tan/3 is only 
allowed to be larger than two as far as we require that the top Yukawa coupling 
constant remains perturbative up to the GUT scale. 

We first show our numerical results for b —>■ s fi + n~ and discuss the electron 
case later. In Fig. ^| and Fig. f|, the branching ratio and the forward-backward 
asymmetry in the SM are shown as functions of the lepton pair invariant mass. In 
the calculation of the b — > s£ + £~ branching ratio we have used m c /mf, = 0.31, 
| V^, / At j = 1.01 and B(6 — > c e V) = 0.104 in Eq. |(3.6)| . We also show similar curves 
for the minimal SUGRA model with a particular set of parameters that tan j3 = 30, 
m = 369 GeV, M gX = 100 GeV, A x = m and the sign of \i is positive. This 
parameter set is chosen so that C-j(m^) has the same magnitude but the opposite 
sign to the SM value. We show the curves with k — ±1 for both models. As 
can be seen in Fig. |3] and Fig. £|, there are large contributions from the J/ip and 
ip' resonances. Since we are interested only in the short distance contribution, we 
consider the following two regions; the low s region, 4mj < s < {mj/^ — 5) 2 , and 
the high s region, {m^> + 5) 2 < s < ml, where 5 is introduced to cut the resonance 
regions and we take 5 = 100 MeV here. We can see that the sizable interference 
between the long and short distance contributions even in these low and high s 
regions. At the asymmetric B factory experiments, however, it may be possible 
to determine the phase of k by measuring the lepton invariant spectrum near the 
resonance regions. Therefore in the following, we consider the branching ratio and 
asymmetry integrated in the above two kinematical region with a choice of k = ±1. 
These are defined as 



B iow(hi g h) = / d§ ( 45 ) 

Jlow(high) 

ds ( / d(cos6) d ^ - f d(cos6) d ^ 
low(high) \J d{cosu)ds J d{cosu)ds 

4°B (hish) = 7^ - r- ( 4 - 6 ) 

^ [ M m d " B f m - d2B 

ds / d{cosu) — + / a(cos( 



low(high) \J d(cos6)ds J d(cos6)ds 



\0 



Notice that for general phase of k the branching ratio takes the value between the 
k = ±1 cases. 
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In Fig. ^, we show the correlation between the branching ratios of the 6 — » s 7 
and b — » s£ + £~ in the above two regions for tan/3 = 30. Fig. |5|(a) and (b) show 
the branching ratio of b — > s£ + £~ in the low s region for k = ±1, and Fig. |5](c) 
and (d) correspond to the high s region. As already mentioned in connection with 
Fig. |2], only (^(m;,) can receive sizable SUSY contributions. It is therefore clear from 
Eq. |(4.1)| and Eq. |(3.6)| that the values of two branching ratios lie on a parabola when 



we neglect SUSY contribution to Cg(mb) and Cio(rrib). This is seen in Fig. [|. If 
we take the experimental constraint on the branching ratio of b — > s 7 into account, 
two separate regions are allowed. One corresponds to the case when the sign of 
CV(^fc) is the same as that of the SM, and the other corresponds to the case with 
the opposite sign. We can see that the branching ratio of b — > s£ + £~ is enhanced 
about 50% in the latter case. Although the branching ratio of b — > s£ + £~ in the 
low s region changes ±15% depending on the sign of k, we can distinguish the sign 
of C-j{mi,) from the branching ratio integrated in this region. On the other hand 
in the high s region, the branching ratio of b — > s£ + £~ depends on the sign of k 
significantly. 

We also show the correlations between the branching ratio of the 6 — » 57 
and the forward-backward asymmetry of the b — > s£ + £~ in Fig. |[ Four figures 
correspond to the case k = ±1 and low/high regions. We can see that the asymmetry 
is also useful to distinguish the sign of CV(m;,). 

We vary the renormalization point Q from to 2m\, in order to study the 
renormalization point dependences, which are also shown in Fig. [| and Fig. |6|. We 
see that the tendency that the branching ratio change along the parabola. This 
is because the change of the renormalization point Q mainly affects Cj(Q). This 
means that we can make a prediction of the branching ratio of b — > s £ + £~ without 
much ambiguity as far as we use the experimental value of the b —>■ 57 branching 
ratio^. Fig. [7| shows that the branching ratio of b — > s£ + £~ in the low s region 
as a function of the chargino mass and the light stop mass for k — 1 and tan/? 
= 30 taking account of the b —>■ sj constraint. The points where the branching 
ratio of the b — > s£ + £~ is enhanced about 50 % compared to the SM correspond 
to the case that the Cf{mb) has the opposite sign to the SM. It is interesting to 

*It is important, however, to reduce the ambiguity of the renormalization point in order to put 
constraints on SUSY parameter space from the b — > s 7 branching ratio. 
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see such parameters correspond to relatively light SUSY particles (m x - < 130 GeV, 
mi < 250 GeV) but beyond the reach of LEP II. We have also analyzed the case of 
small tan/?, for example tan/?= 3. In this case the CV(mft) cannot change its sign 
as shown in Fig. ^ thus the branching ratios change within ±5% after taking into 
account the 5 — > 37 constraint. 

We also calculated the branching ratio and the asymmetry for the b — > s e + e~ 
process. The only difference from the b — > s fi + fi~ case is that the lower limit of 
the lepton invariant mass becomes smaller. Since the CV(m&) term gives dominant 
contribution in the region near the kinematical lower limit, the branching ratio 
and asymmetry integrated in the low s region change from those for b — > s /i + //~ . 
Compared to Fig. || (a), for example, the b — > se + e~ branching ratio is enhanced 
by ~5% at B(6 -> 37) = 1.0 x 10" 4 and ~30% at B(6 ^ s 7 ) = 4.2 x 10" 4 for 
the SM branch and by ~5% and ~20% respectively for the branch with opposite 
sign of CV(mi,). It is worth while noting that we can distinguish the sign of CV(ra&) 
by looking at the low s region in the b —>■ s e + e~ mode just as in the b —>■ s fi + n~ 
case. On the other hand the branching ratio and asymmetry integrated in the high 
s region do not change noticeably from the b —>■ s [i + case. 

Let us now compare our results with those in Ref. |J. As explained in Sec. |T|, 
we have included the one-loop correction term Vj oop in the Higgs potential to find 
appropriate parameter sets. This correction, however, mainly affects the mass of 
the lightest neutral Higgs boson, which does not directly contribute to the FCNC 
processes. Consequently the effect of this improvement is rather smalf| The most 
important difference comes from the long distance contributions of the cc resonances. 
In Ref. 0] b — > s£ + £~ branching ratio is calculated with the short distance contri- 
butions only, omitting the J / if) and if)' resonance regions from the integration range 
of the lepton pair invariant mass. We show that even if the resonance regions are 
avoided, the interference effect between the short distance contributions and the 
tail of the resonances gives ~ ±15% ambiguity to the b — ► s£ + £~ branching ratio 
since the detail of the b-s-J/ip and b-s-tp' couplings, which are parametrized by 
the phase of k in our present analysis, is theoretically unknown. We see that this 
ambiguity is larger than the short distance effects from the SUSY contributions un- 

§ The effect on the fightest Higgs mass wili be important in finding allowed SUSY parameter 
regions when the experimental bound of the lightest Higgs mass is raised. 
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less the CV changes its sign. Thus, it is difficult to extract information about the 
SUSY parameters from the branching ratio and the forward-backward asymmetry 
of b — > s £ + £~ without knowledge of the long distance contributions. This ambiguity 
will be reduced experimentally by the measurement of the behavior of the lepton 
pair invariant mass spectrum around the resonances. This may be achieved before 
the branching ratio in the off-resonance regions is measured since a large number of 
events is expected near the resonance regions. 



V b — > svv decay 

In this section, we present the numerical result of the branching ratio of b — > s v V 
in the minimal SUGRA model. For the calculation of this branching ratio, we need 
to introduce a new operator to the effective Hamiltonian Eq. |(3.1) 



On 



16vr 2 



{sLal^b La ) 



7sM)- 



(5.1) 



i=e,fi,T 



The corresponding Wilson coefficient Cxi is given in [Appendix C| . Note that C\\ does 
not receive QCD correction in LLA. In addition to the SM contribution, there are the 
Z penguin and the box diagrams due to the charged Higgs boson and SUSY particles. 
Since the effect of tau lepton mass in the loop diagram is small, C\\ is calculated with 
3 massless charged leptons. Important difference from the b — > s £ + £r process is that 
no photon penguin diagram can contribute to the b — > s v V process. Thus SUSY 
contributions to C\\ are similar to those to C w , and no large SUSY contribution is 
induced. The branching ratio of b — > s v V is written as 



E B(6 



S Ui Va 



3 B(6 



c e v) 



V, 



cb 



\Cn 



phi 



(5.2) 



In Fig. | we show the scatter plot of the b — > s v V branching ratio and the 
chargino mass (the light stop mass). In this calculation we have taken into account 
all constraints (l)-(7) in Sec. |IV|. We see that the branching ratio does not exceed 
the SM value and the change is at most 10%. This result does not depend much on 
the value of tan /3. 
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VI Conclusions 



In this paper, we have extensively examined the b — > sii branching ratio in the 
minimal SUGRA model. By scanning the three-dimensional space of the soft SUSY 
breaking parameters mo, M g x and Ax for various choices of tan/?, we have found 
that a parameter region, where the Wilson coefficient CV(m&) receives a large SUSY 
contribution, is still allowed under the LEP 1.5 constraints provided that tan/3 is 
large. On the other hand, the SUSY contributions to the coefficients Cg(rrib) and 
Cio(mb) are much smaller than the SM contributions in the whole allowed parameter 
space. Consequently, there is a strong correlation between the predicted values of 
the branching ratios of b — ► sj and b — ► s£ + £~ . Applying the measured bound 
of the b — > 57 branching ratio, we have shown that the predicted values of the 
b — > s£ + £~ branching ratio are separated in two branches for a large tan/3: one 
corresponds to the region where the sign of the b — > sj amplitude is the same as 
that in the SM and the other corresponds to the opposite sign. In the latter case, 
the b — > s£ + £~ branching ratio becomes at most ~ 50% larger compared to the 
SM value. The forward-backward asymmetry is also significantly different from the 
SM in the same parameter region. Since m x - > 100 GeV is allowed for such a 
parameter region, it is possible to observe a large enhancement in the branching 
ratio of b — > s £ + £~ even if no SUSY particle is found at LEP II. We also calculated 
the branching ratio of b — > svv process and found that it is reduced at most 10% 
from the SM value. 

There are several theoretical ambiguities in the calculation of b — > s£ + £~ 
branching ratio, such as the cc resonance effect, the renormalization point depen- 
dence, the strange quark mass and the higher order corrections in the heavy quark 
expansion. Among them we studied the renormalization point dependence and the 
resonance effects in some detail. The strange quark mass correction, which is of order 
rn^/ml, is estimated to be less important especially in the low s region. The higher 
order corrections in the heavy quark expansion is also expected to be small \TE\. 

The renormalization point dependence gives about 30% ambiguity to deter- 
mine the magnitude of Cj{my/) from the measured value of b — > sj branching 
ratio. This ambiguity, however, does not affect the correlation between the branch- 
ing ratios of b — > s 7 and b — ► s £ + £~ much. We can make a rather definite prediction 
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on the value of the b — > s£ + £~ branching ratio with use of the measured 6-^57 
branching ratio. 

The cc resonance effect turns out to be important. To deal with this effect, 
we have introduced a phenomenological parameter k and have presented our results 
for k, — il since the phase of k is not known theoretically. We have pointed out 
that there are sizable ambiguities due to this effect for both low and high s regions. 
This ambiguity will be reduced experimentally if the lepton invariant mass spectrum 
near the resonances will be measured in some detail. We have also found that in 
the low s region the change of the b — > s£ + £~ branching ratio due to the sign of 
C-jirrib) is larger than the ambiguity induced from the phase of k. This fact enables 
us to distinguish the sign of C 7 (m b ) without the knowledge of the phase of k, by 
measuring the b — > s £ + £~ branching ratio integrated in the low s region. 

Acknowledgments 

The authors would like to thank K. Hikasa for carefully reading the manuscript 
and giving useful comments and S. Uno for discussion on the possibility of the 
experimental measurement on the lepton invariant mass spectrum near the resonance 
regions. The works of T. G., Y. O. and M. T. were supported in part by the Grant- 
in-Aid for Scientific Research from the Ministry of Education, Science and Culture 
of Japan. 



16 



Appendix A One-loop Correction to the Higgs Po- 
tential 



The explicit form of Vf oop in Eq. |(2.5)| is given by 



Vi 



loop 
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where Q denotes the renormalization point. The field-dependent masses are given 
by 
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For simplicity we neglect D-term contributions to the scalar masses. 
The radiative electroweak symmetry breaking condition is 



I dV 
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'dV_ 
dh% 



= 0. 



(A.6) 



where the bracket denotes the value at h\ = v cos (3 and h® = v sin f3. We obtain 
the following equations for the SUSY Higgs mass parameter /i and the soft SUSY 
breaking parameter B from |(A.6)| with the explicit form of the Higgs potential |(2.5 " 
and [(XT 1 



v 2 = 



(, 2 + ^ 2 )(tn 2 /3-l) + A?) " ( ^ + A '> tan2 P 
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Appendix B The Effective Hamiltonian for 6 

The effective Hamiltonian for the b ^ s££ is given by 
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v/2 ^ 

in which the operator basis is chosen to be 



01 = (sL Q 7M^a)(cL/37 M CL/3), 

02 = (sL a 7A/3)(cL/37 M C La ), 

03 = (sLal^La) £ (^YqLp), 

q=u,d,s,c,b 

04 = (sLa7A/3) S (VL^qLa), 

q=u,d,s,c,b 
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Appendix C Wilson coefficients at the electroweak 

scale 

In this appendix, we give explicit forms of each contribution to Wilson coefficients 
at the electroweak scale 



Appendix C.l C-j(mw) 



C 7 (m w ) 



cf 



cy + cr +cf +cf + cf, 



-XtXtw 



1 



g/i (a?tw) + h{x t w) 



cot 2 (5\ -fi{x th ) + f 2 {x t h) \ + { T,h{ x th) + U(x th ) 



;\x t h 



Y X Wui(^CLYa2 
q=1 1=1 



XaUjI 3 V Xa u I' 



(C.I) 
(C.2) 

(C.3) 
(C.4) 



+ (T 



,m- 



19 



.2 6 



y y l nob 



(C.5) 
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where Xy = rrqjrn?j and m 8 is the mass of the particle i. fi(x) are the one loop 
functions, which are given in |Appendix E| . The matrix ^cl(r) represents the coupling 
constant of chargino-(up-type-)squark-(down-type-)quark, r^ L( - R ^ represents that 
of neutralino-(down-type-)squark-(down- type-) quark and represents that of 

gluino-(down-type-)squark-(down- type-) quark, which are found in |Appendix D . 
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Appendix C.3 Cg(mw) 
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The matrix ^cur) represents the coupling constant of chargino-sneutrino-lepton 



and r^- L /m that of neutralino-slepton-lepton, which are found in |Appendix D| . Note 
that index % represents the generation of the final lepton and is not summed here. 

Appendix C.4 Cio(mw) 
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Appendix D Feynman rules 



In this appendix, we give our notations. The mass matrix of chargino is given by 

C = -(W~ hT) ( r , ^ a +H.c. (D.l) 

V 1 J \ V2m w cos (3 (j, )\ht) y J 

Diagonalizing this mass matrix, mass eigenstates of chargino Xi(i = 1,2) are given 
by 
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The mass matrix of neutralino is given by 
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where sw = sin 9w, cw = cos 9w, sp = sin (5 and cp = cos (5. Diagonalizing this mass 
matrix, mass eigenstates of neutralino are given by 
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The relevant interaction Lagrangian for the b — > 
follows. 

• chargino-quark(lepton)-squark(slepton) interaction 

£ = ~9Xa 



process is written as 
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The mixing matrices r^j,^) j ^cl(r) <^rid are given by 
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where the matrices Uu, Ue and U v are the unitary matrices which diagonalize 
the up-type squark mass matrix, the slepton mass matrix and the sneutrino 
mass matrix respectively. V is the CKM matrix. Note that we neglect small 
contribution proportional to the Yukawa couplings of the lepton. 



neutralino-quark(lepton)-squark(slepton) interaction 
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The mixing matrices r^ L ^, T £ NL ^ and T V NL are given by 
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where the matrix Ud is the unitary matrix which diagonalizes the down-type 
squark mass matrix. 



d ja d* Ifh (D.17) 



• gluino-quark-squark interaction 

c = -^v / 2(n a ^[(r^)ip L +(r^)iPi i ; 

where the mixing matrices are given by 

iXLYi = (Unfr, (D.18) 
(r^)j = -{Un){ + ". (D.19) 

Appendix E One-loop functions 

These are the one loop functions which appear in calculating the penguin or box 
diagrams. 
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Appendix F The QCD correction 



With Cj(mjy) as the initial condition, we obtain the solution of RGE in the LLA 
approximation |8[,|l4|. 
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where i] 
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+C 2 (m w ) (-0.0185?T - 8994 - 0.071477- 12 / 23 - 0.0380?T - 4230 - 0.0057?? 01456 
-0.42867/ 6 / 23 - 0.6494r ? - 4086 + 2.29 96t? 14 / 23 - 1.0880r/ 6 / 23 ) , 
C 8 (Q) = C^mwW^ (F.8) 

+C 2 (m w ) (-0.0571?T - 8994 + 0.0873r/-°- 4230 + 0.0209?/ 
-0.9135^°- 4086 + 0.8623r/ 14 / 23 ' 



,0.1456 



C 9 (Q) = C 9 (m w ) + , 7r . C 2 (m w ) f-0.1875 + 0.1648r] 1 - - 8994 + 0.242477 1 - 12 / 23 (F.9) 

a s {m w ) \ 

+0.1384r ? 1 - a423 ° - 0.00 737? 1+a1456 - 0.3941r / 1+6 / 23 + 0.0433r ? 1+a4086 ) , 

Cio(Q) = C 10 (m w ), 

C U (Q) = C u (m w ), (F.10) 



a a (m w ) 
M<2) ■ 



Appendix G The kinematical functions 



In this appendix we show the explicit form of kinematical functions in Eq. |(3.6)| and 
Eq. |(3.10)| . Here s, ih s and means s/m 2 , m s /mb and mnjm^. 



w 



(s) = ^/(s-(l + m s ) 2 )(s-(l-m s y), (G.l) 



ai (s } m s ,m e ) = (l + ^H (-25 2 + s(l + m 2 ) + (l-m 2 ) 2 ), (G.2) 
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a 2 (s,m s ,m e ) = (-2s 2 + s(l + m 2 s ) + (1 - m 2 ) 2 ) 



(G.3) 



2mi 



(4s 2 -5(1 +m 2 )s + (1 - 



m 



2\2 



a 3 (s,m s ,m^) 



1 + 



2m? 



(G.4) 



x (-(1 + m 2 )s 2 - (1 + 14m 2 + m 4 s )s + 2(1 + m 2 )(l - m 2 ) 2 ) , 



2m 



a 4 (s,m s ,m,) = 1 + ((1 - m 2 ) 2 - (1 + m 2 )s) . 



(G.5) 
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Figure Captions 



FIG. 1: Feynman diagrams in the SM. (a) b — > 57, (b) b — > sg, (c) the penguin 
diagram for 6 — > s£ + £~, (d) box diagram for b — > s£ + £ _ . The SUSY contributions 
to these diagrams are obtained by replacing the internal lines with SUSY particles. 



FIG. 2: C 7 (m b ), Cg(m b ) and C w (m b ) in the SUGRA model normalized to that of 
in the SM (a) for tan/3 = 3 and (b) for tan/3 = 30. 



FIG. 3: B(6 -> s£ + £-) in the SM and in the minimal SUGRA model for « = ±1. 
The SUSY parameters are fixed with tan/3 = 30, m = 369 GeV, M gX = 100 GeV, 
A x = mo, where C 7 (m b ) becomes the opposite sign to the SM. 



FIG. 4: A F B(b -> s£ + £") in the SM and in the minimal SUGRA model for « = ±1. 
The SUSY parameters are fixed with tan/3 = 30, m = 369 GeV, M g x = 100 GeV, 
A x = m , where C 7 (m b ) becomes the opposite sign to the SM. 
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FIG. 5: A correlation between B(6 — > S7) and B(6 — > s£ + £ _ ) in the minimal 
SUGRA model (a) in the low s region with k — +1, (b) in the low s region with 
/c = — 1, (c) in the high s region with k — +1 and (d) in the high s region with 
k — — 1 for tan/5 = 30. Two vertical dashed line represent the experimental 6—^57 
constraint. Circles, squares and triangles represent how much B( b — > S7) and 
B(6 — > s£ + £~) change when the renormalization point Q is taken to be mb/2,m,b 
and 2m& respectively. 



FIG. 6: A correlation between B(6 — > S7) and *Afb( & - ^ s^ + ^~) in the minimal 
SUGRA model (a) in the low s region with k — +1, (b) in the low s region with 
k — —1, (c) in the high s region with k — +1 and (d) in the high s region with 
k = — 1 for tan/3 = 30. Two vertical dashed line represent the experimental 6^57 
constraint. Circles, squares and triangles represent how much B( b — > 57) and 
-^fb(^ - >■ s^ + ^ _ ) change when the renormalization point Q is taken to be m,f,/2, nib 
and 2mf, respectively. 



FIG. 7: B(6 — > s£ + £ ) in the low s region for tan (3= 30 as a function of (a) the 
light chargino mass and (b) the light stop mass. The solid line shows the value in 
the SM. 



FIG. 8: B(b — > svv) for tan/3 = 30 as a function of (a) the light chargino mass and 
(b) the light stop mass. The solid line shows the value in the SM. 
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Fig. 1 



